home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / DPPCO.z / DPPCO
Text File  |  1996-03-14  |  3KB  |  133 lines

  1.  
  2.  
  3.  
  4. DDDDPPPPPPPPCCCCOOOO((((3333FFFF))))                                                            DDDDPPPPPPPPCCCCOOOO((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      DPPCO   - DPPCO factors a double precision symmetric positive definite
  10.      matrix stored in packed form and estimates the condition of the matrix.
  11.  
  12.      If  RCOND  is not needed, DPPFA is slightly faster.  To solve  A*X = B ,
  13.      follow DPPCO by DPPSL.  To compute  INVERSE(A)*C , follow DPPCO by DPPSL.
  14.      To compute  DETERMINANT(A) , follow DPPCO by DPPDI.  To compute
  15.      INVERSE(A) , follow DPPCO by DPPDI.
  16.  
  17.  
  18. SSSSYYYYNNNNOOOOPPPPSSSSYYYYSSSS
  19.       SUBROUTINE DPPCO(AP,N,RCOND,Z,INFO)
  20.  
  21. DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN
  22.      On Entry
  23.  
  24.      AAAAPPPP DOUBLE PRECISION (N*(N+1)/2)
  25.         the packed form of a symmetric matrix  A .  The
  26.         columns of the upper triangle are stored sequentially
  27.         in a one-dimensional array of length  N*(N+1)/2 .
  28.         See comments below for details.
  29.  
  30.      NNNN INTEGER
  31.         the order of the matrix  A .  On Return
  32.  
  33.      AAAAPPPP an upper triangular matrix  R , stored in packed
  34.         form, so that  A = TRANS(R)*R .
  35.         If  INFO .NE. 0 , the factorization is not complete.
  36.  
  37.      RRRRCCCCOOOONNNNDDDD DOUBLE PRECISION
  38.         an estimate of the reciprocal condition of  A .
  39.         For the system  A*X = B , relative perturbations
  40.         in  A  and  B  of size  EPSILON  may cause
  41.         relative perturbations in  X  of size  EPSILON/RCOND .
  42.         If  RCOND  is so small that the logical expression
  43.         1.0 + RCOND .EQ. 1.0
  44.         is true, then  A  may be singular to working
  45.         precision.  In particular,  RCOND  is zero  if
  46.         exact singularity is detected or the estimate
  47.         underflows.  If INFO .NE. 0 , RCOND is unchanged.
  48.  
  49.      ZZZZ DOUBLE PRECISION(N)
  50.         a work vector whose contents are usually unimportant.
  51.         If  A  is singular to working precision, then  Z  is
  52.         an approximate null vector in the sense that
  53.         NORM(A*Z) = RCOND*NORM(A)*NORM(Z) .
  54.         If  INFO .NE. 0 , Z  is unchanged.
  55.  
  56.      IIIINNNNFFFFOOOO INTEGER
  57.         = 0  for normal return.
  58.         = K  signals an error condition.  The leading minor
  59.         of order  K  is not positive definite.  Packed Storage The following
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. DDDDPPPPPPPPCCCCOOOO((((3333FFFF))))                                                            DDDDPPPPPPPPCCCCOOOO((((3333FFFF))))
  71.  
  72.  
  73.  
  74.      program segment will pack the upper triangle of a symmetric matrix.
  75.         K = 0
  76.         DO 20 J = 1, N
  77.         DO 10 I = 1, J
  78.         K = K + 1
  79.         AP(K) = A(I,J)
  80.         10    CONTINUE
  81.         20 CONTINUE LINPACK.  This version dated 08/14/78 .  Cleve Moler,
  82.      University of New Mexico, Argonne National Lab.  Subroutines and
  83.      Functions LINPACK DPPFA BLAS DAXPY,DDOT,DSCAL,DASUM Fortran
  84.      DABS,DMAX1,DREAL,DSIGN
  85.  
  86.  
  87.  
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.